The Parameterized Post-Priedmann Framework for Theories of Modified Gravity: 

Concepts, Pormahsm and Examples. 
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A unified framework for theories of modified gravity will be an essential tool for interpreting the 
forthcoming deluge of cosmological data. We present such a formalism, the Parameterized Post- 
Friedmann framework (PPF), which parameterizes the cosmological perturbation theory of a wide 
variety of modified gravity models. PPF is able to handle spin-0 degrees of freedom from new scalar, 
vector and tensor fields, meaning that it is not restricted to simple models based solely on cosmo- 
logical scalar fields. A direct correspondence is maintained between the parameterization and the 
underlying space of theories, which allows us to build up a 'dictionary' of modified gravity theories 
and their PPF correspondences. In this paper we describe the construction of the parameterization 
and demonstrate its use through a number of worked examples relevant to the current literature. 
We indicate how the formalism will be implemented numerically, so that the dictionary of modified 
gravity can be pitted against forthcoming observations. 
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I. INTRODUCTION 

Einstein's theory of General Relativity has been re- 
quired to defend its title as the true theory of gravitation 
since its birth. During most of the 20th century modifica- 
tions to General Relativity (GR) were largely an abstract 
venture into the realms of mathematical possibility. To- 
day the issue is a more pressing one that we have been 
forced to entertain by experimental results. 

It is unfortunate, then, that our methods of testing 
GR have become less efficient. Constraining modified 
theories on an individual basis is likely to be an infinite 
process, unless our ingenuity at constructing new theories 
wanes [l| . We need a fast way to test and rule out theories 
if we are to drive their population into decline. 

An analogous situation existed in the 1970s, when the 
question at hand was "Is GR the correct description of 
gravity in the Solar System?" A compelling answer in 
the affirmative was provided by the Parameterized Post- 
Newtonian framework (PPN) [3-11]. In this formalism, 
competing gravitational theories were mapped onto a 
unified parameterization and constrained simultaneously 
using data from laboratory tests, lunar laser ranging and 
early satellite experiments. 

We propose to revive this approach by constructing a 
parameterized framework that can be used to test the 
concept of cosmological modified gravity in a very gen- 
eral, model-independent way (the PPN framework itself 
cannot be applied on cosmological scales, though see Q 
for related ideas). A parametric approach is nothing new; 
there has been a substantial body of work along these 
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lines in recent years However, most approaches 

have considered modifications to the field equations of 
GR that are motivated by simplicity and their relevance 
to a limited number of cases. Model-builders have moved 
beyond simple scalar-field theories, and there is need for 
a parameterization that can handle more sophisticated 
theories. 

In this work we present a new formalism calle d th e 
'Parameterized Post- Friedmann' framework (PPF) [lll| . 
which systematically accounts for the limited number of 
ways in which the Einstein field equations can be mod- 
ified at the linearized level. This means that it encap- 
sulates a very wide variety of theories, without the use 
of approximations or numerical solutions. Table U gives 
a non-exhaustive list of theories that are covered by the 
PPF framework. 

The key feature of our parameterization is that it main- 
tains a direct correspondence between the parameters 
[ll2j | of the formalism and 'known' theories. In this pa- 
per we will refer to a 'known' theory as an established 
model for which field equations can be written down an- 
alytically. These are usually derived directly from a co- 
variant action, though knowledge of the action itself is 
not required for PPF. A known model is represented by 
a point in the space of all possible theories (or a small re- 
gion if the theory contains variable parameters) . In con- 
trast, we will use the description 'unknown' for a point in 
theory-space for which we do not possess the correspond- 
ing action. An unknown theory is characterized purely 
in terms of its PPF parameters. PPF can be used to 
make statements about unknown regions of theory-space 
in addition to the testing of known theories. Such state- 
ments could be of use in guiding model-builders to the 
most relevant regions of theory-space. 

Any parameterization has limits of applicability, and 
PPF is no exception. The (fairly mild) assumptions un- 
derlying our formalism are stated in Table |lTl PPN and 
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Category 


Theory 


References 


Horndeski Theories 


Scalar- Tensor theory 
(incl. Brans-Dicke) 


[18 1 Ql 


f{R) gravity 


[20, 2^ 


f{G) theories 


[22-24] 


Covariant Gahleons 


[25, 26] 


The Fab Four 


[27-30] 


K-inflation 


[3^ 


Generalized G-inflation 


[32, 33] 


Kinetic Gravity Braiding 


[34, 35] 


Quintessence (incl. 
universally coupled models) 




Effective dark fluid 


[38] 


J-JVji-ClluZj V HJlCLLllltk LllCLJl ICO 


Einstein- Aether theory 


[39-42] 


Horava-Lifschitz theory 


[43, 4^ 


> 2 new degrees of freedom 


DGP (4D effective theory) 


[45,46] 


EBl gravity 


[47-5^ 


TeVeS 


[52-5^ 



TABLE I: A non-exhaustive list of theories that are suitable for PPF parameterization. We will not treat all of these explicitly 
in the present paper. Q = — AR^^vR^" + Rfj.upaR^'^'"^ is the Gauss-Bonnet term. 



PPF are highly complementary in their coverage of differ- 
ent accessible gravitational regimes. PPN is restricted to 
weak-field regimes on scales sufficiently small that linear 
perturbation theory about the Minkowski metric is an 
accurate description of the spacetime. Unlike PPN, PPF 
is valid for arbitrary background metrics (such as the 
FRW metric) provided that perturbations to the curva- 
ture scalar remain small. PPF also assumes the validity 
of linear perturbation theory, so it is applicable above the 
scale where matter energy density perturbations cross the 
nonlinear threshold Smikni) ^ 1; note that this boundary 
evolves with redshift. 

Perturbative expansions like PPN and PPF cannot 
be used in the nonlinear, strong-field regime inhabited 
by compact objects. However, this regime can still 
be subjected to parameterized tests of gravity via elec- 
tromagnetic observations [H, [s^ and the Parameter- 
ized Post-Einsteinian framework (PPE) for gravitational 
waveforms [57, 58]. Note that despite the similarity in 
nomenclature, PPE is somewhat different to PPN and 
PPF, being a parameterization of observables rather than 
theories themselves. 

The purpose of this paper is to present the formalism 
that will be used for our future results [11] and demon- 
strate its use through a number of worked examples. We 
would like to politely suggest three strategies for guiding 
busy readers to the most relevant sections: 

i) The casually-interested reader is recommended to as- 
similate the basic concepts and structure of the pa- 
rameterization from mi Al and mi El and glance at 
Table [J to see some example theories covered by this 
formalism. 



ii) A reader with a particular interest in one of the ex- 
ample theories listed in TableUmay wish to addition- 
ally read i|II Bllf f 13] to understand how the mapping 
into the PPF format is performed, and the most rel- 
evant example(s) of mil\ 

iii) A reader concerned with the concept of parameter- 
ized modified gravity in and of itself may also find 
^ff Fl and ijlVI useful for explaining how the approach 
presented here can be concretely implemented (for 
example, in numerical codes). mV\ also discusses the 
connection of PPF to other parameterizations in the 
present literature. 

Our conclusions are summarized in jvl 

We will use the notation k — Mp"^ = 8nG and set 
c — I unless stated otherwise. Our convention for the 
metric signature is (— , +, +, +). Dots will be used to in- 
dicate differentiation with respect to conformal time and 
hatted variables indicate gauge-invariant combinations, 
which are formed by adding appropriate metric fluctua- 
tions to a perturbed quantity (see ijff D[) . Note that this 
means X X- 



II. THE PPF FORMALISM 

A. Basic Principles 

As stated in the introduction, the PPF framework sys- 
tematically accounts for allowable extensions to the Ein- 
stein field equations, whilst remaining agnostic about 
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Assumptions and Restrictions of the PPF Parameterization 

All field equations are second-order or lower in time derivatives (but f{R) gravity is still treatable - see text). 
There exists an FRW solution for the background cosmology.^ 

There exists a frame in which matter components obey their ordinary conservation equations, that is, V^T^ = 0.* 
The field equations of a gravitational theory are gauge form- invariant (see i]IIB|l . 
Nonlinear perturbations are negligible at the lengthscales under consideration. 

If TV non-GR fields are present and N > 2, then TV — 2 relations between the new fields must be specified. If < 2 (the 
majority of cases) then no additional relations are required. 

TABLE II: A summary of the assumptions underlying the PPF formalism. 

^ This potentially poses a problem for Massive Gravity, in which exact, flat FRW solutions do not exist ^6^. However, it may still be 
possible to map the slightly perturbed background solutions of [gH onto our parameterization. 

^ Models that posit a universal coupling between a quintessence field and matter components are treatable, see ^IID3I Models which 
implement non-universal couplings are not. 



their precise form. A particularly important extension 
is to permit the existence of new scalar degrees of free- 
dom (hereafter d.o.f.) that are not present in GR jll3l |. 
We emphasize that scalar d.o.f. are not synonymous with 
scalar fields, but can also arise from the spin-0 perturba- 
tions of vectors and tensors. In addition to these possible 
new d.o.f., there could also be modifications involving the 
d.o.f. already present in GR. Therefore we need to allow 
for perturbations of the spacetime metric to appear in 
the linearized field equations in a non-standard way. 

With these concepts in mind, the construction of the 
PPF parameterization proceeds as follows: first, we add 
to the linearized gravitational field equations of GR new 
terms containing perturbations of all the scalar degrees 
of freedom that are present. Schematically: 

SGf^ = K6Tj^ + 6Ul^ (1) 
where 5U^, = 5U]^^'''-{5gp,) + 5Uf,-°-'-[5^,V,v...) 

In the expression above the tensor of modifications SU^j^u 
has been decomposed into parts containing scalar pertur- 
bations of the metric, scalar perturbations of new d.o.f., 
and matter perturbations. We will see in ijllD 31 that the 
matter term can be eliminated in favour of additional 
contributions to the first two terms of eq.iH]). Possi- 
ble contributions to SV^'"'^' from new vector and tensor 
fields are indicated - schematically for the present - by 
5Ai — {yiV)/a and 2v^ij — —Scjij respectively, where 
7ij is a spatial matrix, is a new vector field and v is 
a spin-0 perturbation of a new tensor field cjpa . It is nec- 
essary to specify the number of new d.o.f. a priori, but 
one can remain completely general as to their physical 
origin. 

Likewise we also need to choose the derivative order 
of the parameterization, that is, the highest number of 
time derivatives that appear in the field equations. We 
will choose this to be two, as Ostrogradski's theorem 
[6^ generically leads to the existence of instabilities in 
higher-derivative theories. Under special circumstances 
these instabilities can be avoided i63||, such as occurs in 
the popular class of /(i?) gravity theories. However, any 



/(i?) theory can be mapped onto a second-order scalar- 
tensor theory via a Legendre transformation [20| . so our 
parameterization is still applicable to /(i?) gravity in this 
form. A reminder of this transformation is given in Ap- 
pendix 1^ 

The coefficient of each perturbation appearing in the 
PPF parameterization can only be a function of zeroth 
order 'background' quantities; for an isotropic and homo- 
geneous universe this means that they are only functions 
of time and scale (through Fourier wavenumber fc). The 
scale-dependence of these coefficient functions is not com- 
pletely arbitrary - it is related to the derivative order of 
theory (see ^IV A|l . We will find that for many known 
theories it has a simple polynomial form. 

Having accounted for all terms permitted in a second- 
order theory, one then imposes gauge form-invariance 
on the parameterization (to be explained momentarily 

in sm. 

One can imagine that the parameterization dictated 
by the above principles could potentially contain large 
numbers of free functions. Fortunately these are not all 
independent; for most cases there exist constraint equa- 
tions that reduce their number considerably. These are 
discussed in j|IIFI below. 

Finally we highlight that as the PPF framework is 
based on the linearized field equations, it is applicable 
to scale regimes that are well above those affected by 
screening mechanisms. Effects such as the chameleon, 
symmetron or Vainshstein mechanisms [63 - [70j are not 
present in the parameterization. 

B. Gauge Form-Invariance 

Any equations of motion that can be derived from 
an action principle must possess the property of form- 
invariance, as we explain here in terms classical mechan- 
ics. 

A canonical transformation is defined to be one that 
takes us from a set of coordinates and their conjugate mo- 
menta {qi,pi}, to a second set {Qi, Pi} in such a way that 
the structure of Hamilton's equations for any dynamical 
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system is preserved. Hamilton's equations are said to be 
form-invariant under canonical transformations. Here, 
form-invariance means that the equations of motion in 
the new coordinate system are simply those obtained by 
replacing {qi,Pi} with {Qi,Pi} in the original system. 
This symmetry is not be confused with general covari- 
ance which, conceptually, enforces that the physics of a 
situation remains independent of coordinate choice. Gen- 
eral covariance allows one to choose a convenient set of 
coordinates such that some perturbation variables are set 
to zero, which will clearly cause the structure of the field 
equations to appear different after the transformation 

One can show |72| that the form-invariance of 
Hamilton's equations carries over to their equiva- 
lent Lagrangian formulation. A gauge transformation 
x^^ ^ x'^ + ^'^ is an example of a canonical transforma- 
tion. Therefore, if a set of gravitational field equations 
are derivable as the Euler-Lagrange equations of an ac- 
tion then they must be form-invariant under gauge trans- 
formations of this kind: they are gauge form-invariant. 

Gauge form-invariance will be realized in our examples 
because any new terms containing the generating vector 
f that arise from a transformation will cancel each other 
by virtue of the background-level field equations, leaving 
the original expression unchanged. See mi El for a demon- 
stration of this. This property will not be 'visible' if one 
has already fixed a gauge. 

To reiterate more concisely: before we choose a gauge, 
the linearized field equations are gauge form-invariant; 
in a specific gauge they are not. Hence enforcing gauge 
form-invariance in our parameterization ensures that it 
does not implicitly correspond to any particular gauge 
choice. We implement this explicitly by combining as 
many terms as possible into gauge-invariant variables, 
and adding a gauge form-invariance-fixing term to guar- 
antee that any remaining gauge- variant pieces cancel via 
the background field equations. Note that this does not 
represent a sleight-of-hand in any sense, as we know that 
the true underlying field equations must be gauge form- 
invariant. 



u!m- Modifications to GR and exotic fluids/dark energy 
are formally indistinguishable at the unperturbed level, 
so X and Y can be regarded as an effective energy density 
and pressure respectively. 



D. Perturbation Variables 

1. Metric Variables 

We write the perturbed line element as (showing scalar 
perturbations only): 



- {1 - 2E)djf - 2{V ^e)dri dx 



dx^ dx^ 



(4) 



where jij is a flat spatial metric. Bardeen introduced 
two gauge-invariant combinations of these perturbations 

^^-^if3-V^iy) + ^H{i> + 2e) (5) 

i' = -E-^{i) + 2e)-^n{i^ + 2e) (6) 

These reduce to the familiar potentials $ and ^E* of the 
conformal Newtonian gauge upon setting e = = - 
though beware different sign conventions. For our pur- 
poses it will be more convenient to use gauge-invariant 
variables of the same derivative order (note that ^ is 
second-order in time derivatives, whilst $ is flrst-order). 
For this reason we define the following combination, in 
which the second-order terms cancel: 

f = i (I. + (7) 

$ and r will be the basic building blocks of the metric 
sector of our parameterization. 



2. New Degrees of Freedom 



C. Background-Level Parameterization 

In eq. ([T]) we introduced a tensor of modifications to the 
field equations of GR. Let the zeroth-order components of 
this tensor be Uq = —X, f/j = YSj, such that the modi- 
fied Friedmann and Raychaudhuri equations are: 

3{H^+K) = Ka'^pM+a^X (2) 

2(H^-n + K'^ = Ka^pM{l+ujM)+a^{X + Y) (3) 

where H is the conformal Hubble factor and dots denote 
derivatives with respect to conformal time. Hereafter we 
specialize to flat cosmologies, setting K — 0. pM denotes 
a sum over standard matter components, including cold 
dark matter, with the effective total equation of state 



We also need to introduce a gauge form-invariant way 
of parameterizing the new scalar degrees of freedom. 
First consider a gauge transformation x^^ ^ x^ -\- ^'^ gen- 
erated by the vector = 1/a (^°, V*V')- The most gen- 
eral way that a dimensionless scalar perturbation x can 
transform under this co-ordinate shift is: 

X ^ X + ^ {Gie + G2i" + Gsij + G4^) (8) 

where Gi — G4 are functions that are flxed in a known 
theory, and we have suppressed their arguments. For ex- 
ample, the perturbation of a standard scalar fleld trans- 
forms as (50 —> 6(f) + /a, so in this instance Gi = (p 
and G2 = G3 = G4 = (we have implicitly normalized (j) 
by a mass scale to reduce clutter whilst keeping x dimen- 
sionless). Table [HT] gives the transformation properties of 
some other common types of scalar degrees of freedom. 
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Perturbation Type 


Symbol 


Gi 


G2 


G3 


G4 




H 





-1 








Metric 


e 


1 





H 


-1 




P 


m 












u 








2 





Fractional energy density 


5 













Velocity potential 


e 


1 











Fractional pressure 


n 


-m{i + uj)cV 











Anisotropic stress 


s 














Scalar field 


5<f) 













Timelike vector: spatial-component* 


V 








n 


-1 



TABLE III: The behaviour of perturbations encountered in this paper under the gauge transformation cc'' — !• a;'' + f , where 
^ = u] — „^,y , — r is the adiabatic sound speed. 

"'"Such that the perturbations to the spatial components of the vector are -written = {S7^V)/a. 



We "want to find a combination of the new d.o.f. and 
metric perturbations such that the gauge transformation 
properties of each term cancel to zero overall. An infinite 
number of such possibilities exists, but we will impose 
the restriction that the resulting gauge-invariant com- 
bination contains no time derivatives. Then our gauge 
invariant variable must have the form: 

X = X + wiS -I- L02P + ^36 + Lo^v (9) 

By using the transformation properties given in eq.® 
and Table Hill and requiring that x has no net transfor- 
mation, one finds: 

wi = G2 a;2 = --^(Gi 4- G4) 

c^3 = G4 LOi = -]-{G^ + HGi)-^{Gi + Gi) (10) 

In the case of a standard scalar field one obtains 
5(1)- ^{P + k'^v). Eqs.® and ^ should be 
thought of a 'template' gauge-invariant variable that can 
be adapted to play the role of the new d.o.f. in a given 
theory. Note that if one uses a perturbation of the space- 
time metric as an input to this prescription the result 
vanishes, since x represents new d.o.f. only. 

3. Matter Perturbations 

A coupling between new d.o.f. and matter can give 
rise to an effective evolving and/or scale-dependent grav- 
itational constant, e.g. in a typical scalar-tensor theory 
Geff = G /4>. This means that when we write the modi- 
fied field equations in the form of eq.((T]) (5C//f will contain 
matter variables. For example, consider a case where the 
modified field equations have the following form: 

^^'^'^ = f(r^ f~ + + (11) 



where once more we have indicated that the modifications 
could arise from new scalar, vector or tensor fields, and 
some arguments have been suppressed. The renormalized 
matter terms are taken into SUf^^: 

(12) 

STf^^ can then be eliminated from the expression above 
using eq.(IT|). Rearranging, we then have SU^^i, expressed 
entirely in terms of metric variables and new d.o.f.: 

SU,. = [1 - /] 5G,, + f [SUX + SUp'^-"^] (13) 

A theory that posits a universal coupling between a 
new d.o.f. and matter species can be transformed to a 
frame where the coupling is absent. These are not nec- 
essarily the standard Einstein and Jordan frames, as the 
field equations do not have to match GR in either frame. 
In the uncoupled frame the the modification tensor U^u 
will be independently conserved, and the PPF param- 
eterization can proceed as normal (we explain why the 
relation V^?7^ = is necessary in §11 Pp . 

In the case of non-universal couplings (s^, |33| it may 
not be possible to find a frame in which the relations 
y^Uj^ = and V^T/f = hold individually. Our param- 
eterization cannot be applied to these cases in its present 
formulation. 



E. Perturbation Framework 

The components of the perturbed Einstein tensor are: 
-a^SG^o = Ea = 2V2$ - 6nkf - m{H^ - n){v + 2e) 
~a^5Gl = V.Se ^2kV+{n^ - ?i)(j> + 2e) 

a^5G\ = Ep^ Qkt + 12-HA:f - 2V^{^ - ^) 

- 6(7^2 _ -^^^ + 3(h^ +UU- (y + 2e) 
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a2,5G} =D}£;e = (14) 
where 

projects out the longitudinal, traceless, spatial part 
of (5G(^ , and we have defined the quantities Ei in order to 
refer to the left-hand side of the Einstein equations with 
ease. We will also use analogous symbols to refer to the 
components of the 5C//f tensor, i.e. 

Up = a^6Ul D]U^ = a'{U]~^5U^5i) (15) 



Having eliminated any matter perturbations from SUj^ 
via the method described in the previous subsection, we 
can expand SUj^ in terms of the gauge-invariant variables 
of eqs.(l5]), (0 and ©• We write down all terms permit- 
ted in a second-order theory. We will see in ^ITFl that the 
equation(s) of motion (hereafter e.o.m.) for the new d.o.f. 
are obtained from the Bianchi identities, and require dif- 
ferentiation of 5Uq and 5Uf\ hence these components can 
only contain first-order time derivatives if the e.o.m.s are 
to remain at second order. 

For simplicity we will consider the case where only one 
new scalar d.o.f. is present, but the extension to multiple 
d.o.f. is straightforward. The PPF parameterization of 
the modified field equations is then as follows: 



Ea = ko^GpmSm + Aak^^ + F^k^f + a^k^x + "i^i + k^M^iv + 2e) 

Ee = na^GpMil + ujm)Om + B^k^ + I^kt + (i^kx + /3ii + k^Me{i> + 2e) 

Ep = 'S ko^GpmUm + Cofc2$ + Cifc$ + Jofc^f + Jifcf + -fok^x + Ti^i + 72^ + k^Mp{v + 2e) 



-r + eox 



k 



X- 



£2_ " 



(16) 
(17) 
(18) 
(19) 



r 



The coefficients ^ ^2 appearing in these expressions 
are not constants - they are functions of the cosmologi- 
cal background, i.e. functions of time and scale. These 
dependencies have been suppressed above for the sake of 
clarity. The factors of k (Fourier wavenumber) accompa- 
nying each term are such that Aq — £2 are all dimension- 
less. A particular known theory will specify exact func- 
tional forms for these coefficients; they can be considered 
as the PPF equivalent of the ten PPN parameters. They 
are the 'slots' one maps a theory of modified gravity onto. 

Ma, Mq and Mp are the gauge form- invariance- fixing 
terms described in mi B[ and are similarly functions of 
background variables. However, the Mi differ from the 
coefficients Aq — £2 in that the former are fixed by the 
zeroth-order field equations, whilst the latter are deter- 
mined by the linearly perturbed field equations. 

Let us demonstrate this explicitly for ea.([T6|. Under 
gauge transformations of the form: 



x^' + C 

ea.([TC| becomes: 

E'A-m{n^ - n) 



= l/a VV) (20) 



Ka pM 

k^MA 



s'M-3n{i + ujM) 



{v' + 2£') 



(21) 



where primed variables indicate those belonging to the 
transformed coordinate system. In order for the form of 



eq. (ITSt to be exactly preserved by the transformation the 
terms containing must cancel, which determines Ma'- 



- -2P 



2(H2-H)-«;aVM(l+c^A/) 



(22) 



where eq. ([3]) has been used in reaching the second equal- 
ity. Analogously one finds: 



Mp = 



2fc3 



(23) 



No gauge form-invariance- fixing term is needed in eq. (|19p 
because all the terms there are individually gauge invari- 
ant, including the matter shear perturbation Sa/ . 

Note that the Mi are unlikely to be useful discrimina- 
tors between theories because their evolution with red- 
shift should be broadly similar in theories that reproduce 
a ACDM-like background. It is in the coefficients Aq — £2 
that the differences between theories will be manifested. 



F. Constraint Equations 



The parameterization laid out in eqs.(|16 |) - (|19|) contains 
twenty- two coefficient functions, — £2. We will see in 
the worked examples of ^Illl that these coefficients are of- 
ten very simply related and not independent. However, 
the purpose of PPF is to obtain non-model-specific con- 
straints from the data - does this mean that we must 
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run a Markov Chain Monte Carlo analysis with twenty- 
two free functions? It would be justifiable to object that 
current and near-future data may not possess sufficient 
constraining power for such a task. 

Fortunately this is not the case for the analyses most 
relevant to current research. Imposing some restrictions 
on the types of theories being considered allows one to 
derive constraint relations between the PPF coefficients, 
thereby immediately eliminating some freedom from the 
parameterization. We will present a set of seven con- 
straints relations here; it is very likely that others exist, 
which we will leave for future investigation (soj . 

This set of seven constraint relations stems from the 
divergenceless nature of the Einstein tensor, V^G(^ = 
(a straightforward consequence of the Bianchi identi- 
ties). We will assume that ordinary matter obeys its 
standard conservation law, V^T^ = 0. Therefore the U- 
tensor must also be divergenceless, V^C//f = 0; the result 
S C^fiUj^) — follows at the perturbative level. Note that 
this statement is not valid for models that involve a non- 
universal coupling between a quintessence field and mat- 
ter species, see i]IID3l 

In an isotropic spacetime the expression S (V^^CZ/f) = 
has two independent components (for v = 0, i) and so 
yields two second-order equations that specify how the 
evolution of x is tied to the metric potentials. How- 
ever, this situation is potentially problematic - how can 
we guarantee that the solutions of these two equations 
agree? In a known theory with a single new field there 
will be one equation of motion for x. How do we rec- 
oncile this fact with the two evolution equations coming 
bomS{V^Uif) = 7 

There are three possible ways that the problem could 
be resolved: 

1. The true e.o.m. is given by the v = component of 
the Bianchi identity. The v = i component reduces 
to a triviality because all of the coefficients of x, $ 
and r in it vanish identically. 

2. Instead the converse is true: the v = i component 
becomes the e.o.m., and the v = Q component re- 
duces to a triviality. 

3. The true e.o.m. corresponds to a combination of 
the u = 0, i components. In this case one set of 
solutions for {x, $, F} would necessarily have to 
solve both equations. 

We will categorize theories as 'type 1', 'type 2' etc. ac- 
cording to which of the three possibilities above occurs. 
All of the single-field theories we tackle in this paper are 
type 1 theories, except for the special case of ^111 Gl as dis- 
cussed below. This is what one might intuitively expect 
- we often think of the e.o.m. of a field as an expression 
of conservation of its energy, and it is the v = Q compo- 
nent of V^r^ = that reduces to the energy conservation 
statement of a classical fluid in Minkowski space. 

One can also construct an argument as to why type 1 
theories are the most natural occurrence by deriving the 



conservation law for U^i, at the level of the action. In this 
derivation the statement V^?7^ = is obtained by con- 
sidering translation along a timelike gauge vector, which 
corresponds to evolving the new fields along a worldline. 
The statement V^J7f = arises from translation along 
the spacelike Killing vectors. As this corresponds to a 
symmetry of the FRW spacetime, it should not lead us 
to dynamical equations. See Appendix [B] for more de- 
tails. 

For a type 1 theory, then, the e.o.m. in terms of the 
parameterization of eqs. ljlGHl^l) is: 



[ai -I- Hfe72] X - 



"0 + A + — + 'Hk{ai + 7i) 



kx 



-^+Hfc(ao+7o)+/3o 



Ao+HfcCi -3-^(X + y) 



k 



k^ ^ + [Fo + HkJi] kt 



y +Hfe(Fo-t- Jo) 



fc^F = 



(24) 



where Hk = H/k. Meanwhile the seven coefficients of 

{<!>, F, F, X, Xi x} ill the second Bianchi component 
= i) must all vanish identically, leading to the follow- 
ing set of constraint equations: 



/3; 



11 + 
3 



re2 =0 



2H/3i)--(7i-2ei) 



^ + 2Hfe/?o-i(7o 



Bo 



- 2eo) = 



^ + 2nkBo-l{Co-2Do)=0 
k 3 

/o - i ( Ji - 2Ki) = 



(25) 
(26) 

(27) 
(28) 
(29) 
(30) 



2HkIo 



i(Jo - 2Ko) 



1 



a'^{X + Y) = (31) 



k " " 3^ " nk 

Thus if one is prepared to focus on type 1 theories - which 
covers a large part of current investigations into modi- 
fied gravity - then seven free functions are immediately 
removed from the PPF parameterization. 

Theories with more than one new d.o.f. cannot be clas- 
sified as type 1/2/3. For a theory with precisely two new 
d.o.f. it is necessary to use both components of the con- 
servation law for U^i, to obtain the required two evolution 
equations. However, if one of the evolution equations can 
be suitably inverted it may be possible to eliminate one of 
the two extra fields; this will modify the PPF coefficients 
of the original system. The new system can then be clas- 
sified as type 1/2/3 if either/neither Bianchi component 
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reduces to a triviality when the new PPF coefficients are 
used. An example of this is given in mil Gl 

For theories with more than two new d.o.f. - for ex- 
ample DGP ( mil C[) and Eddington-Born-Infcld gravity 
( mil Dl) - extra relations between the new fields must be 
provided to close the system of equations. This will pre- 
vent us from being able to constrain regions of unknown 
theory space with more than two new fields. Fortunately, 
based on the relative scarcity of such theories in the lit- 
erature - and simplicity arguments - this type of theory 
is not of primary interest at present. 



always write the above action in a Brans-Dicke-like form, 
with /(</)) = (j) and = uj{(j})/(l) [la]. We wiU work 

with the action in this form. 

As explained in §11 Al and Appendix [Aj a Legendre 
transformation maps f{R) gravity into to a scalar-tensor 
theory, with the following equivalences: 



df{R) 
dR 



= fE 



= 



y(0) = -[i?0-/(i?)] 



(33) 



III. WORKED EXAMPLES 

In this section we demonstrate how a variety of 
commonly-discussed theories of modified gravity can be 
mapped onto our parameterization. We will not dis- 
cuss the motivation or phenomenology of each theory at 
length; our intention is to begin compiling a 'dictionary' 
of theories and their translation into PPF format. 

We will begin with a simple example involving cos- 
mological scalar fields; we progress to more complicated 
cases which demonstrate how theories involving timelike 
vector fields, Lorentz violation and brane scenarios can 
be encapsulated by the PPF parameterization. We then 
treat Horndeski theory, itself a powerful parameteriza- 
tion that subsumes a large portion of theory space. We 
conclude with the example of GR supplemented by an 
exotic fluid. Whilst conceptually simple, this final exam- 
ple has the unusual property that it can be transformed 
from a theory of two new fields to a single-field theory of 
either type 1 or type 2. 



A. Scalar-Tensor Theory and f(R) Gravity 



Hence the expressions below can be straightforwardly 
adapted for use with f{R) models. In fact both theo- 
ries can be obtained as special cases of the Horndeski 
Lagrangian (see ^III F|) . but due to their prevalence we 
will treat them separately here. 

The background effective energy density and effective 
pressure for scalar-tensor theory are (see eqs.© and ([3])): 



a'^X =3-^2(1 - 0) + -uji(t>)-^ - 3n<p + a^Vict)) (34) 
2 (p 

2 



(35) 



The equation of motion for the scalar field is (where mat- 
ter terms have been eliminated using the gravitational 
field equations): 



w(0) 



1 




+ a- 



,dV{^) 



(36) 



A general scalar-tensor theory has an action of the fol- 
lowing form: 



SsT 




-gd^x [ f{4>) R + i^(0)V^0V^0 - 2V{cj)) ] 

(32) 



where ■0" ^re matter fields. For convenience (in this 
section only) we have renormalized the scalar field 
(j) — >■ (j)/Mp such that it is dimensionless. 

In fact one of the functions /((/>) and K{(j)) is redun- 
dant. Through a redefinition of the scalar field one can 



The gauge-invariant combination for perturbations of 
the scalar field dictated by equs.® and ((TU]) is: 



X 



(/? + fcV) 



(37) 



To implement the PPF parameterization, we take the 
linearly perturbed gravitational field equations of scalar- 
tensor theory and regroup terms to form the gauge- 
invariant variables $, F and x- This puts the equations 
into the format of eg s. (| 1611 191) . from which the PPF coef- 
ficient functions can be read off: 



An 



-2(1 

1 
Hk 



4>+ (j) i LJ- 



dv{4>) 



u 



u 

H 



9 



Co = 2(1 



30(3) 30 



30 

mi 



LO (j) 



n 2W 



d(t) 



30 
k7P 



Fn = - 



H{1 



kn<f> 



30 



Jo = — -(1 - 0)(3H + 'iW - k^) - 6-^ - 7:77 2-H - — + -f- 



m 



kU kU 



n 



k^ 



70 



1 dcj(0) 02 

2 d0 "0 

w(0)^-H 



1 dcj(0) 02 

2 d0 

H + w(0)^ 



1 



1 r^^'^' 



3 00 









= 2(1 - 0) - 


Jl 


= 6(1 - 0) - 




= 


ai 





30 
77 



3-H 



72 = 3 

ei = 62 = 



(38) 



One can verify by direct substitution that these coeffi- 
cients obey the constraints of eqs. ([25] -l3i p . Hence scalar- 
tensor theories and f{R) gravity are type 1 theories ac- 
cording to the classification of W Fl 

For computations of observables in f{R) gravity it is 
useful to define a parameter Q, which has the rough inter- 
pretation as the ratio of the Compton wavelength of the 
'scalaron' field fa to the wavelength of a Fourier mode 



Q = 3k 



2 Jrr 



Ac 
A 



(39) 



Making use of the quasistatic approximation (see ^IVBI) , 
the correspondence between our parameterization and Q 
is given by: 



Q 



Do - eo 



l + i^o 

Bo 



Bo 

/3o 



/3o 



This result is most easily obtained by mapping both PPF 
and /(i?) gravity onto the common {/i(a, k), 7(0, k)} pa- 
rameterization discussed in i jlVBI 



B. Einstein- Aether Theory 

Let us now give an example where the new d.o.f. do not 
arise from a scalar field. Modern Einstein- Aether theory 
was introduced by Jacobson & Mattingly [39|] (t hough i n 
fact the earliest incarnation dates back to Dirac |75l - l78| ) 
as a minimalistic model for a theory that violates Lorentz 
invariance. It achieves this by introducing a dynamical 
unit vector field, dubbed the 'aether', that picks out a 
preferred spacetime frame. The aether must be time-like 
at the unperturbed level in order to preserve invariance 
under spatial rotations, and this requir emen t is enforced 
by a Lagrange multiplier in the action jll4l |: 

Sea = ^ j d^^y^ [R + M^K, + \{A" + 1)] 

+SM[r,9f..] (41) 

where M has the dimensions of mass and the kinetic term 
of the aether field is specified by the constants 



(40) Kf, = ci5"^57- + C2<5?<5^ + C3<5^<5^ 



(42) 



For convenience we define a = ci 4- 3c2 + C3. The 
constraint equation from the Lagrange multiplier is 
gi^^A^^A" = -1, which leads to A = (1/a, 0,0,0) in an 
FRW background. 
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In our parameterization of the zeroth-order field equa- 
tions the effective energy density and pressure of the 
aether are: 



1 



a{2H + H' 



(43) 



The aether experiences perturbations around its 
zeroth-order direction. We consider only the spin-0 per- 
turbations here, and write the perturbed aether as 



A = A(°) +SA=-{1 + E, ^V) 



(44) 



where the perturbed constraint from the Lagrange mul- 
tiplier has enforced SA'^ = 'E./a. 

The potential V specifying the perturbation of the 
spatial aether components represents a new degree of 
freedom. Following the prescription of mi D 2i it has 
the transformation properties G4 = —1, G3 — T-L and 
Gi = G2 = 0, so its gauge-invariant partner is: 



V 



V 



= V V -V 

H 2 



(45) 



V has the dimension of a length, so a factor of k is 
included when defining the dimensionless perturbation 
appearing in the parameterization of eqs. (|16ffT9|) . i.e. 
X=kV. 

With the preparations completed, we can now put the 
linearized field equations of Einstein- Aether theory into 
the standard format of our parameterization. The PPF 
coefficients can then be read off: 



^ = c, 1-^ 



Do 



-(ci +C3) 2- 



So = 
Di 



k 

— (Ci +C2 +C3) 



kH V 



n 



(ci + C3) 



k 



n 



Ko=0 

ao = Ukia - ci) 

= (Cl +C2+ C3) 

70 = 2Hka 
72 = 

ei = -(ci + C3) 



Ia=a 
Ji — 3a 

ai = Cl 
/3i =0 



71 = " 

£2=0 



-2Hfc(ci+C3) 



(46) 



These coefficients satisfy the constraints of eas. (|25ll3"T]) . 
hence Einstein- Aether is a type 1 theory. 



C. DGP 

Dvali-Gabadadze-Porrati gravity (DGP) [45] is un- 
doubtedly the most commonly-discussed member of the 
braneworld class of modified gravity theories. It is 
also one of the most extensively tested: recent results 
[79} have placed tight constraints on the normal (non- 
accelerating) branch of the theory in addition to the 
existing constraints on the self-accelerating branch [46| . 
Although increasingly disfavoured, we will present its 
PPF correspondence here as a representative example of 
braneworld scenarios [s^l • 

The theory considers our four-dimensional spacetime 
to be embedded in a five-dimensional bulk, with standard 
matter fields confined to the 4D brane. The ratio of the 
effective gravitational constants in the bulk and brane 
defines a crossover scale, Tc = K5/(2k4) (constrained to 
be of order the horizon scale by supernovae data) , above 
which gravitational forces are sensitive to the additional 
dimension. Below the crossover scale the theory can be 
treated as effectively four-dimensional. However, this is 
not to say that DGP reduces to GR below the crossover 
scale ; the theory allows a new scalar d.o.f. to propagate 
[ll5j | , and enters a scalar-tensor- like regime before strong 
coupling of the scalar switches off departures from GR at 
the Vainshtein scale [t^] ■ 

Neglecting brane tension, the full action of the theory 



1 

2^5 



Sdgp = — / d'^iV^^^ [(5)7?-2A5 

1 



(47) 



+ / d^x 



2k4 



(4)i?, + /:Af(7^^5,.,) 



Cm is the Lagrangian of the brane-confined matter fields 
and '^^•'i? is the Ricci scalar of the 5D spacetime metric 
^^^dtiu- The effective energy density, pressure and equa- 
tion of state that modify the Friedmann and Raychaud- 
huri equations are (see eqs.Q and 



X 
Y 



3e 



-n 



arc 



arcH 
Y 1 



[n + 2%^) 



n 



X 



3 \ 



(48) 
(49) 

(50) 



e = +1 corresponds to the self- accelerating solution 
branch, which is asymptotically de Sitter at late times 
but suffers from ghost pathologies 321 . £ = — 1 corre- 
sponds to the normal solution branch, on which a cosmo- 
logical constant or dark energy is still required to achieve 
accelerated expansion. 

The new fields in DGP arise from the projection of the 
electric part of the bulk Weyl tensor onto the brane. It is 
common to treat the components of the projected Weyl 
tensor £"^1^ as an effective 'Weyl fluid' with a radiation- 
like equation of state, ujw = 1/3. The subdominance of 
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the Weyl terms at late times means that the Weyl fluid is 
usually defined at the perturbative level, and its zeroth- 
order components are neglected. The effective Weyl fluid 
perturbations are: 

q'^Eq = K4a^pwSw (51) 



There are three new fields corresponding to the energy 
density perturbation, velocity potential and anisotropic 
stress of the Weyl fluid, so the two components of the 
four-dimensional [/-conservation law do not provide us 
with sufficient equations of motion to solve the system. 
By considering perturbations of the full five-dimensional 
spacetime one can relate the three Weyl perturbations to 
a single d.o.f., the master variable [s^. However, solution 
of the e.o.m. for the master variable requires gradients 
perpendicular to the brane that are not encapsulated by 
the four-dimensional formalism of this paper. Hence we 
must stick with three new fields, and thus DGP is not 
subject to the classification of i jHFl 

The gauge-invariant Weyl fluid variables are: 

5^ = 5w+^-^^{P + k^^) 



(52) 



The metric PPF coefficients for DGP are: 
3 3 Ko^ pwi^ + uJw)^ 



A, 







Bo 



Co 



Do 



Fo 



lo 



Jo 



Ji 



3 na^pwi^ 



-ujw) 



2kn 



4 + 3 CUE 



3 Ka^piy (1 -I- uJw) 



(2 + 3u;e) 



1 



rlXHk 
3 



(fc2 + iu^ - m) 



rl{X- 
k 



-3Y) 
3 



'Hrl{X 

'rlX 
3 



3y) 



rlXkU 
9 



/c^ + 3-H2(4-|-3lj£;) + 3h) 



Ko = 



k 



HrliX 




sr) 



(53) 



To incorporate three new fields, in principle the template 
of eqs. 



(|16m9p needs to be extended to include two more 
sets of terms similar to the x's, each with their own co- 
efficient functions. However, it turns out that for DGP 
most of these possible extra terms are not needed, and 
the only non-zero ^-coefficients are: 



na pw 



2Xrl fc2 
3 K,(]?pw[^ 



Ldw) 



Xrl 2fc2 
3 3 K,a?pw 



(2 + iuE) 



Xrl 2fc2 

3 na^pwi^ + ^w) 



rl{X + 3,Y) 



fc2 



(54) 



where the superscripts indicate the degree of freedom to 
which each coefficient belongs, and should not be con- 
fused with spacetime indices. Note that, as per the pre- 
vious subsection, factors of k have been used where nec- 
essary to define dimensionless perturbations of the Weyl 
fluid. 

The authors of [13] reported that a distinctive sig- 
nature of DGP is the unusual scale-dependence of new 
terms in the modified Poisson equation. Specifically, they 
found that DGP led to terms linear in wavenumber k as 
the result of a brane-bending mode. PPF is formulated 
in four dimensions, so off-brane terms such as this that 
explicitly probe the additional dimension are not present 
in our coefficient functions. In i jIV Al we will see that 
the four-dimensional theory shares the /c-dependence of 
a scalar-tensor theory. 



D. Eddington-Born-Infeld Gravity 

Let us give another demonstration of how the PPF 
parameterization can handle theories beyond scalar-field- 
type models, this time by considering a bimetric theory. 
Bimetric theories have been the subject of recent intense 
interest in the context of Massive Gravity (see [IHl for 
a review); we will instead focus on the related theory of 
Eddington-Born-Infeld (EBI) gravity, where the coupling 
between the two metrics is substantially simpler. EBI 
gravity possesses a second graviton, but unlike Massive 
Gravity the associated scalar mode propagates and could 
potentially be a ghost. 

EBI gravity arose from the realisation of [131 that a 
theory coupling GR to a Born-Infeld-type Lagrangian 
could be reformulated as an action for two coupled met- 
rics (a closely-related theory based on a bimetric refor- 
mulation of Eddington's original affine action [s^ was 
explored in (iol - fsll lSTl Issj). The bimetric re- writing of 
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the EBI action is 
1 



^EBI 



g{R-2h) + ^-q{K~2\) 

(55) 



Matter couples to the usual spacetime metric g^^i, and K 
is the curvature of the auxiliary metric q^^. A can be 
considered as a cosmological constant for the auxiliary 
metric and I is a length. 

For the purpose of the PPF framework we require the 
cosmological perturbation theory laid out in [48|. First 
note that, since our conformal time coordinate is defined 
by the standard spacetime metric, the 'scale factors' of 
the temporal and spatial parts of the auxiliary metric 
are not constrained to be equal to oneanother nor the 
cosmological scale factor 0(77). Instead we write the un- 
perturbed auxiliary metric as: 



goo = -Z{-qY qi.j Rirffnj (56) 

The background-level effective energy density and effec- 
tive pressure are: 

i?3 



a X — na pebi 



PZa 



a^Y = ko^Pebi = — 



aRZ 



(57) 
(58) 



The scalar perturbations of the auxiliary metric are writ- 
ten in an analogous manner to those of the spacetime 
metric in eq.Q: 

ds^ = -Z{rjf{l-2k)d7f-2R{riY{Vii)d7]dx 



Rivf 



However, we need to be aware that the auxiliary met- 
ric perturbations transform in a slightly different way 
to those of the spacetime metric, due to the possible 
inequality of Z{ri) and R{r]). We must use the same 
gauge transformation vector as for the spacetime met- 
ric, = l/a(^°, VV)- The auxiliary perturbations then 
transform as (using primes to denote the transformed 
variables): 





1 


[(i 




a 






1 


a^Z^ 


e - 


f - 
a 


i?2 




1 


'6R 


/3 


+ - 
a 


R ^ 




2^ 



(60) 



These transformations reduce to those of the or- 
dinary spacetime metric perturbations for the case 



Z{ri) = R{r]) = a{r]). The four scalar perturbations of the 
auxiliary metric are the new fields of EBI gravity. Using 
the now-familiar algorithm of i lIID| the corresponding 
gauge-invariant combinations are: 



D = i) — V 



(61) 



We put the linearized field equations of the spacetime 
metric into the standard format of eqs. (|16m"9|) . with each 

of the gauge-invariant auxiliary perturbations S, e, /3 and 
V acting as a new x-type field. The PPF coefficients are 
(where the superscripts on the Greek coefficients indicate 
to which field they belong): 



A, = - 
Co = 



nk^PZa \ R 
i?3 k ^Z^ 




kH^ZaU \R^ 
iZRa ( R Z 



Hk^P \ R Z 



- 2n 



(62) 



an 



dx' dx^ (59) pke 



kH^Za 

i?3 
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R-" 



2k^PZa 



kH^Za 
3ZRa 



^0 



ZRa 



fj ZRa 
^° " ~2^ 



(63) 



Factors of k are included where necessary to define di- 
mensionless perturbations to the auxiliary metric. The 
four additional fields present in EBI gravity means it is 
not subject to our type 1/2/3 classification. In addition 
to both components of the [/-conservation law, two fur- 
ther relations are required to evolve the system; these are 
provided by eqs.(33) and (34) of gl. 



E. Horava-Lifschitz Theory 



Horava-Lifschitz gravity [43|,|4J] was proposed as a the- 
ory of quantum gravity which achieves power-counting 
renormalizability at ultra-violet energy scales, whilst po- 
tentially recovering GR at low energies. GR is rendered 
non-renormalizable by the scaling properties of the gravi- 
ton propagator, but this situation can be remedied by 
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introducing higher-order derivative terms in the grav- 
itational action [83]- However, the price paid for this 
improved UV behaviour is the appearance of a ghostly 
(negative energy) degree of freedom. The origin of the 
ghost can be traced to the presence of higher-order time 
derivatives in the action. 

Hoi^ava's theory sidesteps this issue by breaking 
Lorentz invariance and treating time and space un- 
equally. The Horava-Lifschitz action contains non- 
linear spatial curvature terms that introduce greater- 
than-second-order spatial derivatives without introduc- 
ing higher time derivatives. Consequently the full diffeo- 
morphism group of GR is broken, and only the smaller 
group of foliation-preserving diffeomorphisms is main- 
tained. As is expected, the breaking of general covari- 
ance introduces a new scalar graviton mode; there has 
been considerable effort made to determine whether this 
mode propagates or not [90l - [9^ . 

After using the ADM formalism to define a preferred 
time direction, two flavours of Horava-Lifschitz gravity 
exist depending on whether the lapse function is al- 
lowed to be a function of both space and time (the non- 
projectable theory) or time only (the projectable theory) . 
In the non-projectable theory the action picks up extra 
terms that depend on the shift vector via: 



6i = ajniV(77,f) 



(64) 



where N is the lapse function and we avoid the standard 
'oj' notation to prevent confusion with the cosmological 
scale factor. 

We will focus on the non-projectable extension of 
Hoi^ava's original action put forward by Bias, Pujolas and 
Sibiryakov 93] , which has dominated much of the recent 
discussion. The Bias et al. theory is motivated by the re- 
quirement that the new scalar mode possesses a healthy 
quadratic kinetic term, and also lifts the 'detailed bal- 
ance' condition of Hoi^ava's original theory. The Bias et 
al. action is: 



5* 



HL 



Ml 



dt d'^x - 



^N[£K-V{g,,,b,)] 



where Ck 



dtd^x^NCM{N,N„g,j) (65) 



1 

2N 



V{g^j,bi) ^-R- cMV- + Mp^{d2R^ + dsRijR'^ 
+ C2b,AV + c^RV ^b' + ...) 
+ Mp^{diRAR + d^V^Rjk'^'-W^ 
+ C46iA26* + C5Ai?V,6* + ...). (66) 

Ellipses indicate terms not relevant for linear spin-0 per- 
turbations, and Ci and di are constants. The recovery of 
GR in the limit A = 1 is not clear-cut, as the new scalar 
mode becomes strongly coupled at low energies. Reso- 
lution of the issue depends on whether the scalar mode 
propagates, see the references above. 



In a flat FRW spacetime the zeroth-order field equa- 
tions of Horava-Lifschitz gravity are identical to those 
of GR but with a rescaled gravitational constant, 
Goff — 2Go/(3A — 1). In terms of our background-level 
parameterization this becomes: 



2v 2 3(1 -A) 

a X = Ka Pm—^ — = 3(1 — A)—— 



(67) 



3A- 1 

^^ = -^(l-A)(2H+H^) (68) 



Moving on to parameterization of the linearized field 
equations, we need to be aware of the reduced dif- 
feomorphism group of Hoi^ava-Lifschitz gravity. The 
theory is invariant under the gauge transformations 
•q — f}{ii), x"^ — x'^(ri,x), i.e. we have lost the ability to 
perform space-dependent reparameterizations of the time 
coordinate. However, by making use of the Stiickelberg 
trick we can restore full general covariance and remove 
the graviton scalar mode from the metric sector of the 
theory; instead we can recast it as a extra field, thereby 
obtaining something similar to a scalar-tensor theory. 

The Stiickelberg trick is implemented by promoting the 
temporal component of the gauge transformation vector 

— l/a(^'', VV) to a new field possessing the necessary 
transformation properties to ensure gauge invariance of 
the whole equation. Once these transformation proper- 
ties are determined we can use our standard procedure 
of mi D 21 to construct the gauge-invariant variable for 
the new Stiickelberg field. Note that we are really using 
a quasi-Stiickelberg trick, as we are leaving the spatial 
diffeomorphisms unaltered. 

Let the dimensionless Stiickelberg perturbation be 
X — (C" has the dimension of length) . The accompa- 
nying gauge-invariant variable for the Stiickelberg field 



X = X + 



(/? + k^v) 



(69) 



For convenience we also define the following higher- 
derivative spatial operators: 



/i 

/2 



2c3 fc2 

2 

C2 e 



Ml a2 



2C5 fc4 



c + 



Ml a2 



C4 fc^ 



(8^2 + Mz) {8di ~ 34) k^ 



Ml 



Mp a 

The PPF coefficient functions are found to be: 



(70) 
(71) 
(72) 



Ao = 3(X-l) + f2{l-^ 



Bo = -{l-X) 
C^o = 3(l-A) (2--^ ) - A ( 1 
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3{H^ -H) 



F. Horndeski Theory 



Fo = 9Hfc(l-A) + -/2 



lo = 3(1 - A) 



1 



9(1-A)CH2+H)-/cVi 



•Hfc 

Ji = 9(1 - A) 
H 2 

«o = [3(A - 1) + A + A] 
/3o = 1 - A 
70 = Hfc [6(1 



A) - /i + M 

-A) 



i^i = 

/5i =0 

71 = 3(1 - A) - /i 
= -y (73) 



The Hofava-Lifschitz coefficients obey the constraint re- 
lations of eas. (|25ll3"T]) . hence it is classified as a type 1 
theory. 

We see that Hofava-Lifschitz theory possesses a dis- 
tinctive signature: nearly all of its PPF coefficients con- 
tain spatial derivative operators that are greater than 
second order. This feature is a result of the particu- 
lar type of Lorentz violation of Hofava-Lifschitz theory, 
which allows time and space coordinates to behave dif- 
ferently under scaling transformations in the UV. 

It is these kinds of distinctive signatures that could 
be a useful tool in guiding us towards (or eliminating) 
particular regions of theory space when confronted with 
data. In the case of Hofava-Lifschitz gravity, the chal- 
lenge will be the detection of strongly scale-dependent 
components of the PPF functions which are likely to be 
subdominant to the 'normal' scale-free and fc^ terms (see 



The last few years have witnessed a resurgence of inter- 
est in Horndeski theory which had lain largely for- 
gotten since 1974 until it was independently re-derived 
by Deffayet and collaborators [s^. Horndeski theory 
is the most general Lorentz-invariant extension of GR 
in four dimensions that can be constructed using a sin- 
gle additional scalar field, with the restriction that the 
equations of motion must remain second order in time 
derivatives. All theories that fit this description can be 
obtained via special choices of four arbitrary functions of 
the scalar field that appear in the Horndeski Lagrangian. 
A non-exhaustive list of theories that fall under the Horn- 
deski umbrella is; Brans-Dicke and scalar-tensor gravity, 
f{R) gravity (in its scalar-tensor formulation), single- 
field quintessence and K-essence theories, single-field in- 
fiation models, the covariant Galileon, the Fab Four, 
Dirac-Born-Infeld theory. Kinetic Gravity Braiding, ac- 
tions involving derivative couplings between a scalar field 
and the Einstein tensor and f{Q) theories, where Q is the 
Gauss-Bonnet term. 



Inevitably, the price to be paid for this powerful gener- 
ality is that Horndeski theory is cumbersome to calculate 
with. Its action is 1951: 



Shd = j d^x^/^ A -|-£m(V'",3m!^)^ 



(74) 



where 



C2 = K{cj),X) 

£3 = -G3(<^,A')n</) 

£4 = G4(</.,Af)i?-f G4,;t [(□</-)' -(V^V,0)(V^V'' 0)] 



(75) 
(76) 
(77) 



£5 = G5(0,A')G^,(V^V^ 



G5,;t [{n4>f - m<t>) (V^V,</)) {V^'Vcf) + 2(V^V„0) (V"V^0) (V^V^(/.)] (78) 



K and Gi (i=3,4,5) are four functions of a scalar field and its kinetic energy, X = —d^cj)df^4>/2. GR is recovered by 
setting G4 = Mp/2 and G3 = G5 = K = 0. We will also need to define the following quantities [95l |: 



Ft = 2 



Gt 



Ga-X\ -^G 



G4 - 2XG 



4,X 



5,X 



X 



— t- G5 X + G5 



G 



5,X 



G 



5,0 



(79) 
(80) 
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Gs{<l),X) 


04(0, X) 


G5(0,A') 


Scalar-tensor theory 












f{R) gravity 


-^{RfR-fiR)) 





— JR 





the covariant GaUleon 




■^x 


2 




Horndeski's original notation 
(used for the Fab Four) 




— 2(6ki,00 + 8^ti_0 






Kinetic Gravity Braiding 


K{<i>, X) 


G3(0, X) 







Quintessence & phantom fields 


eX - V{<i>) 










K-essence & K-inflation 


K{^,X) 











TABLE IV: This table lists the choices for the four free functions of the Horndeski Lagrangian that reproduce some previously- 
studied theories of modified gravity. In the table above, if a function is left in general terms then the choice is arbitrary. For 
the covariant Galileon, d are dimensionless constants and m is a mass scale. For the fourth line, kf f Ki dX and expressions 
for the Hi that give rise to the Fab Four are presented in [SJ. In the penultimate line, e = +1 for quintessence and e = — 1 for 
phantom scalar fields. 



(5XG5.X + 2X^G5, 



XX ) 



2'HX (36*5,0 + 2XGti,(j>x) 



(81) 



T = a'iXK^x + ^X^K^xx) + 12'H^XG^^x + m^X^G^^xx ~ 2a^{XG^^^ + X'^G^^^x] 



m^Gi + 6 



-H^ {7XG4,x + WX^Gi^xx + iX'^Gi^xxx) ~ (64^0 + bXGi^4,x + 2X^Gi^ 



4>xx) 



[30n^XG5,x + 2m^X^Gz^xx + m^X^G^^xxx] - 6H'A'(6G5^0 + ^XGz,4,x + 2X^G. 



5,4>XX j 



(82) 



where the relation Xdx4> — 4>l'^ has been used. Note that we have used conformal time, contrary to the authors of 
[95} . The gauge- invariant perturbation of the scalar field is: 



Y — (p + k u) 

^ Mp Mp en ^ ^ 



(83) 



Using a tilde to denote division by the square of the reduced Planck mass (ie. Qt = Gr/Mp = kGt), the PPF 
coefficients for Horndeski theory are: 
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So = — na^pM - 2{n^ ^'H) — 



kU 



n 
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•HP 



4-H 



2 1 

H H 



3 



20(^) [Gr-^4 



fc2 



1202 

k^ 

2Gt- 



Gt I 

/fc2 [ 


2n + 
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] /j2-^2 


ko^Pm 


Igt 






QHa^ 


3a2V(/. 






•H2fc2 


7^A;2 



2aQ 
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Gr 
H 



[I-Gt) 




-2k^{l - Gt) + Sna^PM - 6-^(ae) + 6{n^ +H)- 6^(2^^ - 



drj 



J.=6(l-- 



Ko = -^{I^Gt) 



Q!q = Mp 

Mp 



/3o 



70 



71 



2Mp 
Mp 



l,^-j{ea-nGr 

-20 (ae - HGt) - 2GTH(f> + 2ae-H0 



Ki=0 

2Mp 



- Mp 



ko^Pm 



kcj) 



Gt + nGr - 'HTr ) + iMp-^V 



ai 



Pi = 



72 



k4 

2Mp 



Qa - HGt 



'-^(Ba-nGr 



6^T ( H + 2%^ - 2ntj + 6— (ae - HGt) + GaQ \^3n - 24 j - Sku'^pm 

ei = £2 = 



eo = ^ Itr + nGr - n^r 



(84) 



fi and V are derivatives of the zeroth-order field equations 
with, respect to the scalar field (00 and ii components 
respectively) - see [95| for the relevant expressions. 

Table IIII Fl collects some 'settings' for the Horndeski 
Lagangian functions that reproduce theories of current 
interest. The application of the PPF formalism to Horn- 
deski theory immediately brings a large realm of theory 
space within reach of our parameterization. 



G. GR with a Dark Fluid 

Our final example should really be classed as a dark 
energy model rather than a theory of modified gravity, 
though arguably the distinction is not important. We 
present the example of an adiabatic dark fluid character- 
ized by a constant equation of state lud and negligible 
anisotropic stress. We will use this as an example of how 
a theory with two additional fields can be recast as a 
single- field theory (and hence type 1/2/3) in some cases. 
The example also has relevance to theories that can be 
usefully written as an effective fiuid at the level of the 
linearized gravitational field equations, e.g. quintessence 
and its progeny (whilst an effective fluid interpretation is 
possible for all theories, it is not always useful). 

The zeroth-order modiflcations to the fleld equations 



r 



are simply: 

a^X — na^pp) 



na lodPd 



(85) 



The two new flelds are the fractional energy density 
perturbation and the velocity perturbation, Sp) and 9p), 
which appear in the [/-tensor as: 



Ua = ho^pdSd 
Up — Sko^ pdujdSp) 



C/e =0 



(86) 



The fluid velocity is given by — V^Od- For an adia- 
batic, shear-free fluid the conservation and Euler equa- 
tions are: 



5d = -{1 + iop) [k^eo + - k^e 



(87) 



3d - -H(l - iu:D)OD 



1 + ujd 



These correspond to the two components of the conser- 
vation equation V^tZ/f = 0. Eq. (|87| can be used to elimi- 
nate 6d in favour oi 5d from the gravitational field equa- 
tions, resulting in a theory with a single new field. After 
forming the relevant g.i. combination 



5D = 5D + \{l+uj){P + k'v) 



(89) 
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the PPF coefficients can be extracted (with ah those not 
stated being zero); 



Pi 



(90) 



where il^i is the ratio of the energy density of the dark 
fluid to the critical density pc = iTi^ /{no?), as per the 
usual definition. The superscripts indicate that these co- 
efficients 'belong' to the 5d perturbations, and should 
not be confused with spacetime indices. The e.o.m. is 
obtained by using eq. (j87| to replace 0d in eg. (j88| . lead- 
ing to: 



5D+'H5D{l-2,uJD) + k^i^D So + iil + uoD)^) 

^3 ^ ^2 . 



(91) 



Using the PPF coefficients of the modified single-field 
system, we find that this is our first example of a type 
2 theory. That is to say, the coefficients of each term in 
ea.p4| vanish identically. 

Of course we can implement an analogous procedure to 
eliminate 5d and treat 9d as the single extra field. Note 
that this cannot be done for the special case ojd = 0, for 
which we cannot invert eg. ([88)) to find 5d{0d,So) (and 
using eq. (|87| instead would lead to integral expressions 
that do not fit into our parameterization). The relevant 
gauge-invariant quantity in this case is; 

^D = OD~^{P + k^v) (92) 

Repeating similar steps to before we obtain the PPF co- 
efficients for the new theory, where now Ojj is the only 
additional field. Using tildes to indicate that these are 
not the same as eqs. ((90)) . the results are: 



Bo = -mk^D{l+OJD) 







— SHk^D 




Jo 






~k§ 

"o 


^■mlnD[i-iojD)— 


+ ^d) 

UJD 


-fee 












if 




+ W_d) 


if 


^mlnoii + ujo) 





This time the 0-component of the Bianchi identity pro- 
vides the e.o.m.: 



1 



(94) 



(93) 



The constraints of eas. (l25ll3"T|) are all satisfied. Hence in 
this formulation the dark fluid becomes a type 1 theory. 

We have at last found a counter-example to the 
monopoly of type 1 theories. But is this an 'artificial' 
result of moulding a two-field theory into single-field for- 
mat? We discuss the issue a little further in Appendix 

m 



IV. THE PPF COEFFICIENTS 

In the previous section we demonstrated that a wide va- 
riety of modified gravity theories could be put into the 
standard format of the PPF framework, parameterized 
by a set of coefficient functions. This operation would 
have no real benefit if the results were identical for all 
theories. Fortunately we find that different theories leave 
characteristic signatures in the set of coefficient functions 
that in principle could be used to distinguish between 
them. In the first part of this section we consider one 
aspect of these signatures - their scale-dependence - and 
how this should guide the implementation of PPF in nu- 
merical calculations. 

We will not discuss here the required accuracy needed 
to achieve these distinctions between theories or the fea- 
sibility of obtaining it with current or planned experi- 
ments; we leave that to a future work [5^. The most 
exciting tests of modified gravity in the next decade or so 
are likely to come from measurements of the growth rate 
of large-scale structure, probed via weak gravitational 
lensing, redshift space distortions and peculiar velocity 
surveys . The lengthscales at which these observables 
are relevant lie in the 'quasistatic' regime, so we will dis- 
cuss the behaviour of the PPF framework in this limit in 
iVBl 

In the final part of this section we briefiy connect our 
work to a recent parameterization of screening mecha- 
nisms 0, which are designed to suppress the effects of 
modified gravity on small scales. 



A. Signatures and Scale-Dependence 

A theory of modified gravity leaves its signature in 
both the magnitude and scale-dependence of the PPF 
coefficients. Some of its PPF coefficients may be identi- 
cally zero (see the previous section for numerous exam- 
ples), but otherwise their magnitudes will be controlled 
by the fundamental parameters of the theory. These 
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TABLE V: The scale-dependence of the terms appearing in the modified Poisson equation (|95[) for the example theories of 
ijllll Recall that F oc k~^ , see eq.((7|. Note that in some cases the scalar perturbation x includes a factor of k so that it is 
dimensionless, eg. x ~ in Horava-Lifshitz gravity, = fc6 for a fluid velocity potential. EBI and DGP gravity will contain 
multiple copies of terms 3 and 4 (eg. xi = ^w, X2 = kOw for DGP), but both terms have the same scale-dependence. 



could be adjustable parameters that can be fine-tuned 
to fit growth-rate data, or they may already be tightly 
constrained by other considerations (fitting to the back- 
ground expansion rate, the need to pass Solar System 
tests of gravity, etc.). It is difficult to make general state- 
ments about the magnitudes of the PPF coefficients, as 
these are theory-specific. 

Instead we will focus on their scale-dependence. Usu- 
ally the full modified Einstein equations are not needed 
to calculate observable quantities (power spectra); often 
only the Poisson equation, the slip relation ($ — \I>), and 
the evolution equations for 5m and 9m are used. The 
modified Fourier-space Poisson equation in the PPF pa- 
rameterization is (taking ea. ([TC|) +3"H x eg . (fT7)l ) : 

-2fc^$ = Ka^pM^M 

+ iAo + mkBo)k''^ +{Fo + mkio) k^r 

+ {ao + iUkPa) fc'x +(ai + iUkPi) k i (95) 

where Am = 6m + 3H(1 + ijJm)9m is a gauge-invariant 
density perturbation. 

Table IIVI indicates the scale-dependence of the non- 
standard terms of eq. ([M)) in each of the theories of mill 
It is clear that modifications to the Poisson equation have 
the form of a power series in even powers of the Fourier 
wavenumber k. It is worth stressing that this form is ex- 
act, and not a Taylor series expansion of a more compli- 
cated function; an even power series in k is the only pos- 
sibility. In most cases only the scale-free and quadratic 
terms are present, in agreement with the findings of [s^l 
(apart from DGP, as discussed in mil Cp . This is no sur- 
prise - we have focused on theories that are second order 
in time derivatives, and Lorentz invariance then implies 
that they must be second order in spatial derivatives too. 
Odd powers of k are forbidden by parity. 

Hof ava-Liftschitz gravity is a notable exception - its ex- 
plicit breaking of the symmetry between time and space 
coordinates allows higher-order spatial derivatives to en- 
ter the PPF coefficients. Einstein- Aether theory is dis- 
tinguished by its quadratic components in terms 2 and 



4. 

Another distinctive exception is the k~'^ in term 4 of 
the dark fluid scenario {Sd case). This has arisen from 
the replacement of 9o hy Sd using the energy conserva- 
tion equation, which gives (schematically) 9d ~ Sjj/k'^. 
Note that the k~^ term does not correspond to an intrin- 
sically non-local modification of gravity; such apparently 
non-local terms can be generated by the elimination of 
degrees of freedom from a theory. 

Vanishing or purely scale-independent corrections are 
the calling card of EBI gravity. The scale-dependence 
of scalar-tensor theory and /(i?) gravity must necessar- 
ily be contained within that of Horndeski theory, and is 
largely shared by DGP. Again, this is not surprising since 
DGP passes through a scalar-tensor-like regime on scales 
below the crossover scale Tc, before transitioning to the 
Vainshtein regime at short distances. 

Interestingly, all the theories except Hoi^ava-Lifshitz 
gravity are found to have scale-independent slip relations 
(for brevity we have not decomposed the slip relation into 
its individual terms in Table ITVT) . The coefficients Ki and 
£2 are zero in all cases. 

Based on these observations there are two ways to im- 
plement the PPF coefficients in an Einstein-Boltzmann 
solver code such as the CAMB package [97| : 

i) We can construct a sensible parameterization of 
the PPF coefficients by splitting them into purely 
time-dependent functions multiplying a spatial 
dependence. To be fully general (at least in terms of 
the theories treated in this paper) we should expand 
each of the PPF functions appearing in the Poisson 
equation as: 

4 

Ao = J2Mv)k'^"-'^ (96) 

n=0 

and similarly for Bq, Fq etc. One must choose a 
motivated ansatz for the time-dependent function 
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Aoir]), such as a Taylor series in i}\ or {'Ho/'H)~^. 
In reality, a form as general as ea.(P51) is probably 
unnecessary for all of the coefficients. 

ii) Alternatively one may employ a Principal Co mpo - 
nent Analysis (PCA), as advocated in [lol [gj-fTooj . 
The PCA approach expands the functions to be 
constrained in a basis of orthogonal eigenmodes in 
{k, zj-space, using the eigenbasis that a given data 
set is most sensitive to. Informally speaking, PCA 
reveals 'what the data knows', rather than 'what we 
would like to know'. 

Method i) has the advantage of simplicity, and it uti- 
lizes the common features of modified gravity theories 
to jump directly to the most relevant regions of theory 
space. However, by choosing a form for functions such 
as ^0(^7) we are imposing our preconceptions on the be- 
haviour of modified gravity, namely that it must emerge 
at late times (z < 1) in order to reproduce the effects of 
an apparent dark energy. 

Method ii) maintains a greater degree of agnosticism, 
avoiding the use of semi-arbitrary functional forms. How- 
ever, it effectively throws away much of our prior knowl- 
edge about the structure of physically-reasonable theo- 
ries. One might argue that method i) is optimized for 
constraining known regions of the modified gravity land- 
scape, and method ii) is better suited to exploring new 
unknown territory. 



B. The Quasistatic Limit 

On sufficiently subhorizon scales {Hk <C 1) - but above 
the nonlinearity scale - the time derivatives of metric 
perturbations are small relative to their spatial deriva- 
tives, and can be neglected (at least for the purposes of 
N-body simulations). This is known as the quasistatic 
regime [g^. It is believed that time derivatives of the 
additional fields in modified gravity theories can be sim- 
ilarly neglected in this regime: many scalar-field type 
models exhibit damped oscillatory solutions below a spe- 
cial lengthscale, which rapidly decay into insignificance. 
Only the averaged, stationary spatial profile then remains 
of relevance. 

Let us consider type 1 theories, as these are the most 
prevalent. In the quasistatic limit ea. ((24)) . the e.o.m. for 
the new scalar d.o.f., reduces to a relation between the 
perturbations of the extra field and the potential $: 



treatment to the slip relation, neglecting the anisotropic 
stress of matter. The results can be written in the form: 



[do + k(3o] X + 



Ao + kBo 



$ = 



(97) 



Time derivatives of the PPF functions have been main- 
tained because the quasistatic approximation only allows 
us to neglect time derivatives of perturbations. To ob- 
tain the quasistatic limit of the Poisson equation we use 
eq. (|97l) to replace x in eq. (j95|) . and drop the x and T 

terms (recaU that r = l/k{^ + H^), where $ will be 
small and is suppressed by Hfc). We apply a similar 



where 



Ka'^n{a, k) pm^m 




(98) 
(99) 

(100) 
(101) 



We see that on quasistatic scales modifications to the 
Poisson and slip relations can be reparameterized in 
terms of two time- and scale-dependent functions, /i(a, k) 
and j{a,k). This form of parameterization has been 
exp lored in a large number of investigations 
[T2I [3 [l^ , using a number of different ansatzes for the 
scale-dependence of {/i,7}. From the discussion of the 
previous subsection we see that one does not really have 
the freedom to choose just any ansatz here; with a few 
exceptions, the physically-relevant choice is that of scale- 
free terms plus k^ terms, such as used in [7]. 

Deviations of {/i, 7} from their GR values of {1, 1} 
are degenerate between modifications to the metric sector 
and new d.o.f. sector, indicated by the presence of Greek 
and Latin PPF coefficients in both. The advantage of 
maintaining the unnapproximated PPF Poisson equation 
feg. dMI) ') is that these coefficients are kept distinct, which 
is necessary if we hope to detect the signatures of particu- 
lar theories discussed in the previous subsection. Keeping 
terms such as x in the equations may be unfeasible or un- 
necessary for N-body simulations, but it should not pose 
a problem for Einstein-Boltzmann codes. Furthermore, 
lifting the quasistatic approximation means that the PPF 
parameterization is valid all the way up to horizon scales, 
which are relevant for large-scale CMB modes contribut- 
ing to the late-time Integrated Sach s- Wolfe effect and the 
lensing-ISW cross-correlation [lOl| . 



C. Connection to Screening Parameterizations 

Brax and collaborators have introduced a new param- 
eterization [96i . il02j that is optimized for screeni ng m ech- 
anisms, such as the chameleon [G^lGSj. dilaton |l03j and 
symmctron mechanisms (gtI . l68l|. In this parameteriza- 
tion a theory is described by the time-evolving mass of 
the scalar field, m(a), and its coupling to matter in the 
Einstein frame, /3(a), where the Jordan and Einstein met- 
rics are related by: 

and the coupling function is: 
dliiA{(f)) 



(102) 



/3(a) 



(103) 
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The parameterization is constructed explicitly in the 
Einstein frame; this means that when mapped to our 
(Jordan-frame) parameterization it mixes PPF coeffi- 
cients from the metric and extra-field sectors. In the 
quasistatic regime the mapping is: 



2 



2 



l-Do + eo 



ao + fc/3o 
'Ao_+kBo 



2/32 



= 1 



2/32 



(104) 
(105) 



Adding the two expressions above gives an additional 
constraint relation between the PPF coefficients in the 
quasistatic regime. 

If we lift the quasistatic approximation then in princi- 
ple there is sufficient information to express the twenty- 
two PPF functions in terms of /3(a) and 771(0). How- 
ever, the expressions are not particularly elegant and we 
will not p resen t them here on grounds of relevance: re- 
cent work 104| indicates that the scalar fields involved in 
the chameleon, dilaton and symmetron screening mech- 
anisms must have a Compton wavelength < 1 Mpc in 
order to satisfy Solar System and Galactic constraints. 
Deviations from GR are then only expected in the non- 
linear regime, to which PPF does not apply. 



V. CONCLUSIONS 

We have presented a new framework, the Parameter- 
ized Post-Friedmann formalism (PPF), for conducting 
model-independent tests of modified gravity. The con- 
struction of the framework does not rely on ad-hoc mod- 
ifications to GR, but is built by considering the limited 
number of ways that the linearized Einstein equations 
can be extended whilst maintaining the properties of a 
physically-relevant gravitational theory. An exact, ana- 
lytic mapping exists between a theory of modified gravity 
that obeys the assumptions of TablelUand the set of PPF 
coefficient functions. 

There are four key reasons we believe that the PPF 
formalism represents a significant step forward in param- 
eterized treatments of modified gravity: 

1. New degrees of freedom, a common feature of many 
models, are explicitly parameterized for. Until re- 
cently j96.] this has not generally been the case. 
This is not a minor issue; new degrees of freedom 
are at the heart of the phenomenology that might 
render a theory of modified gravity distinguishable 
from ACDM. 

2. The parameterization is not restricted to models 
based solely on scalar fields. As experimental pre- 
cision increases, the continued success of ACDM 
has forced model-builders to look to richer theories 
that additionally involve vector and tensor fields. 



PPF is able to handle the spin-0 perturbations of 
these theories as easily as canonical scalar fields. 

3. Once cast into the format of PPF, a theory of 
modified gravity may leave a distinctive 'calling 
card' in the set of PPF functions. Such signa- 
tures may allow some theories to be ruled out, or 
at least indicate the general characteristics of the- 
ories that are consistent with the data, eg. one ex- 
pects pure scalar field models to share similarities, 
whilst Lorentz- violating theories may share similar 
traits. It seems reasonable to suppose that higher- 
dimensional models may have their own unique fea- 
tures. 

4. The parameterization is valid from horizon scales 
down to the scales at which nonlinearities become 
important. This spans a wide range of observ- 
ables t argeted b y the next generation of experi- 



105l4l08l |. including weak lensing (of both 



ments 

galaxies and the CMB), redshift space distortions, 
peculiar velocity surveys, the ISW effect and associ- 
ated cross-correlations, and possibly galaxy clusters 

m- 

The PPF framework can be applied in two modes de- 
pending on the user's interests. It can be used as a tool 
for tackling observations, a gateway that allows data con- 
straints to be applied to a multitude of 'known' theories 
simultaneously. By 'known' theories we mean existing 
models for which we have the field equations in hand. Al- 
ternatively, one can use the framework as an exploratory 
tool for model-building. As mentioned in 3) above, PPF 
may be able to indicate the physical features that lead 
to tension with the data, thereby guiding theorists to the 
regions of theory space that are likely to prove most fruit- 
ful. We note that for theories involving more than two 
new fields the exploratory mode cannot be used, but the 
observational mode remains fully functional. 

The cost of the increased capabilities of the PPF frame- 
work is the introduction of more 'free' functions than 
other parameterizations (though we have seen in this pa- 
per that the structure of these functions is most definitely 
not arbitrary) . In ijllFI we presented a number of con- 
straints that can be used to reduce this freedom, and we 
believe that there exist further relations not presented 
here [s^. The situation is likely to be analogous to that 
of the Parameterized Post-Newtonian framework, where 
a set of ten new potentials and ten new parameters are 
necessary to capture nearly all possible distortions of the 
GR metric: it turns out that many of these new param- 
eters are needed in only a handful of special cases, and 
are zero the rest of the time. 

We have been ambitious. It is possible that the cur- 
rent and next-generation datasets may not be able to 
tease apart the degeneracies between contributions to the 
modified field equations to the extent required here. In 
this situation the PPF framework still functions as a null 
test for non-GR behaviour, even if it cannot reach its full 
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potential as a discriminator of modified gravity theories. 
Then the most pertinent question to ask is: how much 
better do our experiments need to be to reahze this dis- 
criminatory power? Tackhng this question requires im- 
plementation of the PPF framework in a numerical code 
for computing cosmological observables. The results of 
this implementation - and a possible answer to our ques- 
tion - will be presented in a future work j59|] . 
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Appendix A: Transformation of f{R) gravity into a 
Scalar- Tensor Theory 

Consider the action: 



S 



d'x ^ [/(C) + /'(C)(i? - C)] + SM{r,9^.) 



(Al) 



where primes denote differentiation with respect to the 
new scalar Variation of this action gives the field equa- 
tion for C as: 



/"(C)(i?-C) = o 



(A2) 



Hence C = i? for aU /"(C) ^ 0. Substituting R = ( into 
eq. (IAl[) recovers the standard f{R) action. Therefore 
ea. (IAl[) is an equivalent action for f{R) gravity, with the 
special case of /"(C) = corresponding to the Einstein- 
Hilbert action. 

To recast the equivalent action as a scalar-tensor the- 
ory we make the definitions: 



df{R) 



dR 

Then eq. (|Al|) becomes: 



/ d^xy^[<j>R~2V{cb)]+SM{i^'',9^..) 



(A3) 



(A4) 



which has the form discussed in ijlll Al with uj{(f>) = in 
this case. 

Note that because the metric has not been transformed 
its couplings to the scalar field and matter are unaffected. 



Appendix B: Is the Bianchi 0-Component Always 
the E.O.M.? A Plausibility Argument 



In ^11 Fl we introduced a classification scheme for single- 
field theories, according to whether the spatial or tempo- 
ral part of the perturbed conservation law SCVfiUj^) — 
reduced to a trivial relation (types 1 and 2 respectively), 
or neither component did (type 3). Of the seven examples 
considered in ijllll we found that four were type 1 (Scalar- 
tensor, Einstein- Aether, Hofava-Lifshitz and Horndeski). 
EBI and DGP (the effective 4D theory) were are not 
subject to this classification, having more than two new 
fields. The adiabatic, shear-free dark fluid could be recast 
as either type 1 or type 2. 

It seems, then, that type 1 theories dominate. Why is 
this? We have yet to find an example of a type 3 theory, 
and it could be argued that our example of a type 2 
theory is the result of an artificial re-writing. We present 
here a plausibility argument why all single-field theories 
may be type 1, but do not make any claims that it is a 
fully rigorous or watertight proof. 

Consider the following action: 



S 



R 

2k. 



C 



M 



-Cr 



I —g d X 



(Bl) 



where Lu is the Lagrangian that gives rise to the U- 
tensor that appears in the gravitational field equations, 
eq.(IT]). We assume that any gravitational action can be 
written in the form of eg. dBip by suitable choice of 
e.g. for /(i?) gravity we would have L\j — ^[/(i?) — R\- 
The variation of the last term in eq. ljBip is: 



u 



1 

2k, 



where 

U^"" = - 



~9 Sg^„ 



(B2) 



(B3) 



Now consider an infinitesimal co-ordinate transforma- 
tion, x'^ —?' x^^ + S,^ (x) . The corresponding change in the 
metric is given by Jg^^ — 
this expression into ea. (|B2 
between indices ^ and i^: 



V^^!, - Vi,^^. Substituting 
and utilizing the symmetry 



SSu = - (V^e.) C/^'^y^d'a; (B4) 



(B5) 



The first term of ea. (|B5|) can be converted into a bound- 
ary integral via the usual divergence theorem. This 
boundary term must equal zero, because it is an assump- 
tion of the variational procedure that the small variations 
vanish at spatial infinity. 
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Demanding stationarity of the action, the second term 
of eq. (|B5p yields the expression: 

rV^C/,^ - (B6) 

which is easily recognized as the conservation law for 
U'^'^ contracted with the gauge transformation vector. 
We obtain the two components of the conservation law 
for Ufii, by considering two separate coordinate transfor- 
mations, a temporal one (x*^ — ^ + and a spatial 
one (x* — > a;* + ^*). The first corresponds to transla- 



tion along a worldline, and the second to translation of 
spatial coordinates across a fixed-time hypersurface. A 
vector proportional to is not a Killing vector of the 
FRW spacetime, so it is not surprising that a coordinate 
translation along results in an evolution equation for 
the new degrees of freedom. In contrast, spatial vectors 
proportional to are Killing vectors of the FRW metric - 
so no dynamics can be obtained from a translation along 
Hence the terms of the ly — i equation must vanish 
identically. 



T. Clifton, P. G. Ferreira, A. Padilla, and C. Skordis, 
Physics Reports 513, 1 (2012). [26 
K. S. Thorne and C. M. Will, Astrophys. J. 163, 595 
(1971). 

C. M. Will, Astrophys. J. 163, 611 (1971). [27 
C. M. Will, Living Reviews in Relativity 9 (2006), URL 
http : //www. livingreviews . org/lrr- 2006-3 
H. Noh and J.-c. Hwang, ArXiv e-prints (2012), [28 
1206.5034. 

W. Hu and I. Sawicki, Physical Review D 76, 104043 [29' 

(2007) . 

E. Bertschinger and P. Zukin, Phys. Rev. D 78, 024015 

(2008) , 0801.2431. [30; 
E. V. Linder and R. N. Cahn, Astroparticle Physics 28, 

481 (2007), arXiv:astro-ph/0701317. [31 

C. Skordis, Physical Review D 79, 123527 (2009). 

L. Pogosian, A. Silvestri, K. Koyama, and G.-B. Zhao, 
Phys. Rev. D 81, 104023 (2010), 1002.2382. [32 
S. F. Daniel and E. V. Linder, Phys. Rev. D 82, 103523 
(2010), 1008.0397. 
R. Bean and M. Tangmatitham, Phys. Rev. D 81, [33 
083534 (2010), 1002.4197. 
T. Baker, P. G. Ferreira, C. Skordis, and J. Zuntz, Phys. [34 
Rev. D 84, 124018 (2011), 1107.0491. 
J. N. Dossett, M. Ishak, and J. Moldenhauer, Phys. Rev. 
D 84, 123001 (2011), 1109.4583. [35; 

D. B. Thomas and C. R. Contaldi, Journal of Cosmology 

and Astroparticle Physics 12, 013 (2011), 1107.0727. [36; 
J. Zuntz, T. Baker, P. G. Ferreira, and C. Skordis, 
Journal of Cosmology and Astroparticle Physics 6, 032 [37 
(2012), 1110.3830. 

R. A. Battye and J. A. Pearson, Journal of Cosmology [38 
and Astroparticle Physics 7, 019 (2012), 1203.0398. 
C. Brans and R. H. Dicke, Phys. Rev. 124, 925 (1961). [39; 
G. Esposito-Farese and D. Polarski, arXiv.org gr-qc, 
63504 (2000). [40; 
T. P. Sotiriou and V. Faraoni, Reviews of Modern [41 
Physics 82, 451 (2010). 

A. De Felice and S. Tsujikawa, Living Rev.Rel. 13, 3 [42 
(2010), 1002.4928. 

S. Nojiri, S. D. Odintsov, and M. Sasaki, Phys. Rev. D [43 

71, 123509 (2005), arXiv:hep-th/0504052. [44 

L. Amendola, C. Charmousis, and S. C. Davis, Journal 

of Cosmology and Astroparticle Physics 12, 020 (2006), [45 

arXiv:hep-th/0506137. 

T. Koivisto and D. F. Mota, Physics Letters B 644, 104 [46' 
(2007), arXiv:astro-ph/0606078. 

A. Nicolis, R. Rattazzi, and E. Trincherini, Phys. Rev. [47 



D 79, 064036 (2009), 0811.2197. 

C. Burrage, C. de Rham, and L. Heisenberg, Journal 
of Cosmology and Astroparticle Physics 5, 025 (2011), 
1104.0155. 

C. Charmousis, E. J. Copeland, A. Padilla, and 

P. M. Saffin, Phys. Rev. Lett. 108, 051101 (2012), URL 

http : //link . aps . org/doi/10 . 1103/PhysRevLett . 10 8.051101[ 
C. Charmousis, E. J. Copeland, A. Padilla, and P. M. 
Saffin, Phys. Rev. D 85, 104040 (2012), 1112.4866. 
J. -P. Bruneton, M. Rinaldi, A. Kanfon, A. Hees, 
S. Schlogel, and A. Fiizfa, ArXiv e-prints (2012), 
1203.4446. 

E. J. Copeland, A. Padilla, and P. M. Saffin, ArXiv e- 
prints (2012), 1208.3373. 

C. Armendariz-Picon, T. Damour, and V. Mukhanov, 
Physics Letters B 458, 209 (1999), arXiv:hep- 
th/9904075. 

T. Kobayashi, M. Yamaguchi, and J. Yokoyama, 
Progress of Theoretical Physics 126, 511 (2011), 
1105.5723. 

C. Deffayet, X. Gao, D. A. Steer, and G. Zahariade, 
Phys. Rev. D 84, 064039 (2011), 1103.3260. 
C. Deffayet, O. Pujolas, I. Sawicki, and A. Vikman, 
Journal of Cosmology and Astroparticle Physics 10, 026 
(2010), 1008.0048. 

R. Kimura, T. Kobayashi, and K. Yamamoto, Phys. 

Rev. D 85, 123503 (2012), 1110.3598. 

L. Amendola, Phys. Rev. D 62, 043511 (2000), 

arXiv:astro-ph/9908023. 

V. Pettorino and C. Baccigalupi, Phys. Rev. D 77, 
103003 (2008), 0802.1086. 

W. Hu, Astrophys. J. 506, 485 (1998), arXiv:astro- 
ph/9801234. 

T. Jacobson and D. Mattingly, Physical Review D 70, 
24003 (2004). 

T. Jacobson, ArXiv e-prints (2008), 0801.1547. 

T. G. Zlosnik, P. G. Ferreira, and G. D. Starkman, Phys. 

Rev. D 75, 044017 (2007), arXiv:astro-ph/0607411. 

T. G. Zlosnik, P. G. Ferreira, and G. D. Starkman, Phys. 

Rev. D 77, 084010 (2008), 0711.0520. 

P. Hofava, Phys. Rev. D 79, 084008 (2009), 0901.3775. 

T. P. Sotiriou, Journal of Physics Conference Series 283, 

012034 (2011), 1010.3218. 

G. Dvali, G. Gabadadze, and M. Porrati, Physics Let- 
ters B 485, 208 (2000). 

W. Fang, S. Wang, W. Hu, Z. Haiman, L. Hui, and 
M. May, Phys. Rev. D 78, 103509 (2008), 0808.2208. 
M. Baiiados, Phys. Rev. D 77, 123534 (2008). 



23 



M. Bafiados, P. G. Ferreira, and C. Skordis, Phys. Rev. 
D 79, 063511 (2009). 

M. Bafiados and P. G. Ferreira, Physical Review Letters [82 
105, 011101 (2010), 1006.1769. 

P. P. Avelino and R. Z. Ferreira, ArXiv e-prints (2012), [83 
1205.6676. 

C. Escamilla-Rivera, M. Banados, and P. G. Ferreira, [84 
Phys. Rev. D 85, 087302 (2012), 1204.1691. 
J. D. Bekenstein, Phys. Rev. D 70, 083509 (2004), [85 
arXiv:astro-ph/0403694. 

C. Skordis, Phys. Rev. D 74, 103513 (2006), URL [ 
http : //link . aps . org/doi/10 . 1103/PhysRevD . 74 . 103513 

C. Skordis, Phys. Rev. D 77, 123502 (2008), URL 

^http : //link . aps . org/doi/10 . 1103/ PhysRevD . 77 . 123502 [87 
T. Johannsen and D. Psaltis, Astrophys. J. 716, 187 
(2010), 1003.3415. [ 
T. Johannsen and D. Psaltis, ArXiv e-prints (2012), 
1202.6069. [89 
N. Yunes and F. Pretorius, Phys. Rev. D 80, 122003 
(2009), 0909.3328. ' [90; 

N. Cornish, L. Sampson, N. Yunes, and F. Pretorius, 
Phys. Rev. D 84, 062003 (2011), 1105.2088. [91 
T. Baker, P. G. Ferreira, and C. Skordis, in prep. 
G. D'Amico, C. de Rham, S. Dubovsky, G. Gabadadze, [92 

D. Pirtskhalava, and A. J. Tolley, Phys. Rev. D 84, 
124046 (2011), 1108.5231. [93 
M. S. Volkov, ArXiv e-prints (2012), 1205.5713. 
M. Ostrogradski, Mem. Acad. St Petersbourg 6 (1850). [94 
R. P. Woodard, Lecture Notes in Physics, vol. 720 of 
Lecture Notes in Physics (Springer, Berhn, Heidelberg, [95 
2007). 

J. Khoury and A. Weltman, Physical Review Letters [96 
93, 171104 (2004). 
J. Khoury and A. Weltman, Physical Review D 69, [97 
44026 (2004). 

J. Khoury, arXiv.org 1011, 5909 (2010). [98; 
K. Hinterbichler and J. Khoury, Physical Review Letters 
104, 231301 (2010). [99; 
K. Hinterbichler, J. Khoury, A. Levy, and A. Matas, 
Physical Review D 84, 103521 (2011). 

P. Brax, C. van de Bruck, A.-C. Davis, B. Li, [lOO; 

B. Schmauch, and D. J. Shaw, Physical Review D 84, 
123524 (2011). [101 
A. L Vainshtein, Phys. Lett. B 39, 393 (1972). 

K. A. Mahk and D. R. Matravers, ArXiv e-prints (2012), [102 
1206.1478. 

L. Hand and J. Finch, Analytical Mechanics (Cambridge [103 
University Press, Cambridge, England, 1998). 
J. M. Bardeen, Physical Review D (1980). [104 
L. Pogosian and A. Silvestri, Phys. Rev. D 77, 023503 
(2008), 0709.0296. [105; 
P. A. M. Dirac, Nature 168, pp.906 (1951). 
P. A. M. Dirac, Proc. R. Soc. Lond. A 209, pp.291 [106' 

(1951) . 

P. A. M. Dirac, Proc. R. Soc. Lond. A 212, pp.330 [107 

(1952) . 

P. A. M. Dirac, Proc. R. Soc. Lond. A 223, pp.438 

(1954). [108; 

A. Raccanelli, D. Bertacca, D. Pietrobon, F. Schmidt, 

L. Samushia, N. Bartolo, O. Dore, S. Matarrese, and [109 

W. J. Percival, ArXiv e-prints (2012), 1207.0500. 

R. Maartens and K. Koyama, Living Reviews in Rela- [110 

tivity 13, 5 (2010), 1004.3962. 

C. Charmousis, R. Gregory, N. Kaloper, and A. Padilla, [111 



Journal of High Energy Physics 10, 066 (2006), 
arXiv:hep-th/0604086. 

D. Gorbunov, K. Koyama, and S. Sibiryakov, Phys. Rev. 

D 73, 044016 (2006), arXiv:hep-th/0512097. 

S. Mukohyama, Physical Review D (Particles 62, 84015 

(2000). 

R. Wagoner and R. D. Blandford, Monthly Notices of 
the Royal Astronomical Society (2008). 
K. Hinterbichler, Reviews of Modern Physics 84, 671 
(2012), 1105.3735. 

A. Eddington, The Mathematical Theory of General 
Relativity (Cambridge University Press, Cambridge, 
England, 1924). 

J. Casanellas, P. Pani, I. Lopes, and V. Cardoso, Astro- 
phys. J. 745, 15 (2012), 1109.0249. 
Y.-H. Sham, L.-M. Lin, and P. T. Leung, ArXiv e-prints 
(2012), 1208.1314. 

I. Kimpton and A. Padilla, Journal of High Energy 
Physics 7, 14 (2010), 1003.5666. 

C. Charmousis, G. Niz, A. Padilla, and P. M. Saffin, 
Journal of High Energy Physics 08, 070 (2009). 

D. Bias, O. Pujolas, and S. Sibiryakov, Journal of High 
Energy Physics 10, 029 (2009). 

X. Gao, Y. Wang, R. H. Brandenberger, and A. Riotto, 
Physical Review D 81, 83508 (2010). 

D. Bias, O. Pujolas, and S. Sibiryakov, Physical Review 
Letters 104, 181302 (2010). 

G. W. Horndeski, International Journal of Theoretical 
Physics 10, 363 (1974). 

A. De Felice, T. Kobayashi, and S. Tsujikawa, arXiv 
gr-qc (2011). 

P. Brax, A.-C. Davis, B. Li, and H. A. Winther, Phys. 

Rev. D 86, 044015 (2012), 1203.4812. 

A. Lewis, A. Challinor, and A. Lasenby, Astrophys. J. 

538, 473 (2000), arXiv:astro-ph/9911177. 

G.-B. Zhao, L. Pogosian, A. Silvestri, and J. Zylberberg, 

Physical Review Letters 103, 241301 (2009), 0905.1326. 

A. Hojjati, G.-B. Zhao, L. Pogosian, A. Silvestri, 

R. Crittenden, and K. Koyama, Phys. Rev. D 85, 

043508 (2012), 1111.3960. 

G.-B. Zhao, R. G. Crittenden, L. Pogosian, and 
X. Zhang, ArXiv e-prints (2012), 1207.3804. 

E. Di Valentino, A. Melchiorri, V. Salvatelli, and A. Sil- 
vestri, ArXiv e-prints (2012), 1204.5352. 

P. Brax and P. Valageas, ArXiv e-prints (2012), 
1205.6583. 

P. Brax, C. van de Bruck, A.-C. Davis, and D. Shaw, 

Phys. Rev. D 82, 063519 (2010), 1005.3735. 

J. Wang, L. Hui, and J. Khoury, ArXiv e-prints (2012), 

1208.4612. 

Dark Energy Surve y (DES), 

|http : // www ■ darkenergy survey ■ or g [ 

Large Synoptic Survey Telescop (LSST), 
http : //www. Isst . org/lsst' 

A. Refregier, A. Amara, T. D. Kitching, A. Rassat, 
R. Scaramella, J. Weller, and f. t. Euclid Imaging Con- 
sortium, ArXiv e-prints (2010), 1001.0061. 
Square Kilometre Array (SKA), 

http : //www . skatelesco pe . org[ 

M. Tegmark, Phys. Rev. D 66, 103507 (2002), 
arXiv:astro-ph/0101354. 

K. Koyama, Phys. Rev. D 72, 123511 (2005), URL 
http : / /link . aps . org /doi/10 . 1103/PhysRevD . 72 . 12 35111 

We note that this is not the first time the phrase 'Pa- 



24 



rameterized Post-Friedmann' has been used [6l. ll09l |. We 
have chosen to recycle the name here (with kind permis- 
sion from W. Hu) because we beheve an analogy with 
PPN is a concise and broadly accurate way of describing 
our formalism. However, we must stress that the analogy 
with PPN should not be taken to extremes - hopefully 
the differences between PPN and PPF are made clear 
in this paper. 

[112] We will see in the next section that in an expanding 
universe one is forced to use time-dependent functions 
as 'parameters', rather than constant numbers. 

[113] This is the most common method that mainstream mod- 
ified gravity theories use to evade Lovelock's theorem. 
Lovelock's theorem states that the only gravitational 
theory with second-order equations of motion and no 
new degrees of freedom (beyond those of GR) that can 



be derived from a covariant action in four dimensions is 
that of GR with a cosmological constant term. 
Ifowever, some theories instead evade Lovelock's theo- 
rem by invoking more than four dimensions, nonlocality, 
or higher-derivative field equations. 

[114] In this paper we treat only the linear Einstein- Aether 
theory. In the generalized Einstein-Aether theory K, is 
replaced by an arbitrary function T{JC). This will lead 
to more complicated PPF coefficients but is not intrin- 
sically problematic. 

[115] The new scalar d.o.f. can arise from either a brane- 
bending mode (if the brane tension is negative), the 
helicity-0 mode of the graviton (if the brane tension is 
positive), or a su perposition of the two (for zero brane 
tension) P.[Tlol]. 



